Emittance Growth from the Thermalization of Space- Charge Nonuniformities* 
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Abstract 

Beams injected into a linear focusing channel typically 
have some degree of space-charge nonuniformity. In gen- 
eral, injected particle distributions with systematic charge 
nonuniformities are not equilibria of the focusing channel 
and launch a broad spectrum of collective modes. These 
modes can phase-mix and have nonlinear wave-wave inter- 
actions which, at high space-charge intensities, results in a 
relaxation to a more thermal-like distribution characterized 
by a uniform density profile. This thermalization can trans- 
fer self-field energy from the initial space-charge nonuni- 
formity to the local particle temperature, thereby increasing 
beam phase space area (emittance growth). In this paper, 
we employ a simple kinetic model of a continuous focus- 
ing channel and build on previous work that applied system 
energy and charge conservation[[j] g] to quantify emittance 
growth associated with the collective thermalization of an 
initial azimuthally symmetric, rms matched beam with a 
radial density profile that is hollowed or peaked. This emit- 
tance growth is shown to be surprisingly modest even for 
high beam intensities with significant radial structure in the 
initial density profile. 

1 INTRODUCTION 

Experiments with high-current, heavy-ion injectors have 
observed significant space-charge nonuniformities emerg- 
ing from the source. Sharp density peaks on the radial 
edge of beam have been measured, but the local incoher- 
ent thermal spread of particle velocities (i.e., the particle 
temperature) across the beam is anticipated to be fairly uni- 
form since the beam is emitted from a constant temperature 
surface. When such a distribution is injected into a linear 
transport channel, it will be far from an equilibrium condi- 
tion (i.e., particles out of local radial force balance), and a 
broad spectrum of collective modes will be launched. 

The spatial average particle temperature of a heavy ion 
beam emerging from an injector is typically measured as 
several times what one would infer from the source ther- 
mal temperature (~ O.leV) and subsequent beam envelope 
compressions, with f x ~ 20eV where f x ~ [e 2 x /{2R 2 )]£ b . 
On the other hand, the radial change in potential energy 
from beam center to edge is qA<fi ~ 2.25keV for a beam 
with line-charge density A ~ 0.25/iC7m(A</> ~ A/(47reo))- 
If even a small fraction of such space-charge energy is ther- 
malized during collective relaxation, large temperature and 
emittance increases can result. 

In this paper, we employ conservation constraints to bet- 
ter estimate emittance increases from collective thermal- 
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ization of normal mode perturbations resulting from ini- 
tial space-charge nonuniformities characteristic of intense 
beam injectors. Past studies have employed analogous 
techniques to estimate emittance increases resulting from 
the thermalization of initial rms mismatches in the beam 
envelope and space-charge nonuniformities associated with 
combining multiple beams and other processes 111 0, Bfl. 

2 THEORETICAL MODEL 

We analyze an infinitely long, unbunched (d/dz = 0) non- 
relativistic beam composed of a single species of particles 
of mass m and charge q propagating with constant axial 
kinetic energy £ b . Continuous radial focusing is provided 
by an external force that is proportional to the transverse 
coordinate x, i.e., F rxt = —2£hk 2 Rn x., where krm = const 



bK^O 3 



is the betatron wavenumber of particle oscillations in the 
applied focusing field. For simplicity, we neglect particle 
collisions and correlation effects, self-magnetic fields, and 
employ an electrostatic model and describe the transverse 
evolution of the beam as a function of axial propagation 
distance s in terms of a single-particle distribution function 
/ that is a function of s, and the transverse position x and 
angle x' = dx/ds of a single particle. This evolution is 
described by the Vlasov equation[Q], 
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(1) 



fc 2 x 2 /2 + (q/2£b)(p is the single- 



where H = x' 2 /2 

particle Hamiltonian and the self-field potential <f> satisfies 

the Poisson equation (CGS units here and henceforth) 



V 2 = -Airq / d 2 x' f 



(2) 



subject to the boundary condition <f>(r = r p ) = at the 
conducting pipe radius r = |x| = r p = const. 

If no particles are lost in the beam evolution, the Vlasov- 
Poisson system possesses global constraints corresponding 
to the conservation of system charge (A) and scaled energy 
(U) per unit axial length, 



A 



U 



d 'rx Id 'fx' f = const, 
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Here, W = Jd 2 x |V0| 2 /(8tt) is the self-field en- 
ergy of the beam per unit axial length and (£) = 
(Jd 2 x Jd 2 x' £ / )/(Jd 2 x Jd 2 x' / ) is a transverse statis- 
tical average of £ over the beam distribution /. Note that U 
includes both particle kinetic energy and the field energy of 
the applied and self-fields. These conservation laws follow 
directly from Eqs. (|l|)-(||) and provide powerful constraints 
on the nonlinear evolution of the system. 



Moment descriptions of the beam provide a simplified 
understanding of beam transport. For an azimuthally sym- 
metric beam {d/d6 = 0), a statistical measure of the beam 
edge radius R = 2{x 2 ) 1 / 2 is employed. Note that R is the 
edge radius of a beam with uniformly distributed space- 
charge. Any axisymmetric solution to the Vlasov-Possion 
system will be consistent with the rms envelope equationfnj 



ds 2 + ho* R R3 - ° • 



(4) 



Here, Q = qX/Eb = const is the self-field perveance and 



4[(x 2 )(x K 
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] ' is an edge measure of the 



rms x-emittance of the beam and is a statistical measure of 
the beam area in x-x' phase-space (i.e., beam quality). For 
general distributions, e x is not constant and evolves accord- 
ing to the full Vlasov-Poisson system. 

3 NONUNIFORM DENSITY PROFILE 

We examine an beam with an azimuthally symmetric radial 
density profile n = Jd 2 x' f given by 



n(r) 




< r < r b , 
r b <r < r p . 



(5) 



Here, r b is the physical edge-radius of the beam, uq = 
n(r = 0) is the on-axis (r = 0) beam density, and h and p 
are "hollowing" [0 < h < oo, h = n(r = rb)/n(r = 0), 
p > 0] and radial steepening parameters associated with 
the density nonuniformity. This density profile is illustrated 
in Fig. [I] for the steepening index p — 2 and hollowing fac- 
tors h — 1 (uniform), h = 1/2 (hollowed), and h = 2 
(peaked). The hollowing parameter h has range < h < 1 
for an on-axis hollowed beam and < 1/h < 1 for an 
on-axis peaked beam. The limit h — ► 1 corresponds to a 
uniform density beam and h, 1/h — ► correspond to hol- 
lowed and peaked beams with the density approaching zero 
on-axis and at the beam edge (r = rj,), respectively. For 
large steepening index p ^> 1, the density gradient will be 
significant only near the radial edge of the beam (r ~ r&), 
and the density is uniform for h — 1 regardless of p. 
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Figure 1: Uniform, hollowed, and peaked density profiles. 

The beam line-charge density (A) and rms edge-radius 
(R) are related to the parameters in Eq. (|j) by 

X = I a x n = irqnor b 



(p + 2)h 



R = 2{x 2 ) 1 ' 2 



(p + 4)(pft + 2) 



n 



(6) 



Using these expressions, the Poisson equation (M) can be 
solved for the potential <j> corresponding to the density pro- 
file (@) and used to calculate the self-field energy W as 



W = X 



4(p- 
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It is convenient to define an average phase advance pa- 
rameter a for the density profile (g|) in terms of an enve- 
lope matched (R' = = R"), rms equivalent beam with 
uniform density (h = 1) and the same perveance (Q) and 
emittance (e^) as the (possibly mismatched) beam with a 
nonuniform density profile (h ^ 1). Denoting the phase 
advance per unit axial length of transverse particle oscilla- 
tions in the matched equivalent beam in the presence and 
absence of space-charge by a and ctq, we adapt a normal- 



ized space-charge parameter cr/cro = \/kp — Q/R 2 /fc/3o- 

The limits <j/<Jo — ► and cr/cro — > 1 correspond to a cold, 
space-charge dominated beam and a warm, kinetic domi- 
nated beam, respectively. Note that this measure applies 
only in an equivalent beam sense. In general, distributions 
/ consistent with the density profile (||) will not be equi- 
libria (d/ds 7^ 0) of the transport channel and will evolve 
leaving a ill defined. 

4 EMITTANCE GROWTH 

We consider an initial beam distribution / with a density 
profile given by Eq. (||) and an arbitrary "momentum" dis- 
tribution in x'. Such an initial distribution is not, in general, 
an equilibrium of the focusing channel and a spectrum of 
collective modes will be launched (depending on the full 
initial phase-space structure of /). These modes will phase- 
mix, have nonlinear wave-wave interactions, etc., driving 
relaxation processes that have been observed in PIC simu- 
lations to cause the beam space-charge distribution to be- 
come more uniform for the case of high beam intensities. 
The conservation constraints (g|) are employed to connect 
the parameters of an initial (subscript i), nonuniform den- 
sity beam with /i ^ with those of a final (subscript /), 
azimuthally symmetric and rms envelope matched beam 
(R' f = = R" f ) with uniform density (h = 1). 

Employing Eqs. (Ph-M), conservation of charge (Xi — 
Xf = X) and system energy (Ui = Uf) can be combined 
into an single equation of constraint expressible as 

(i?//i?,) 2 -l , P (l-h)[4 + p+(3+p)h] 



l-(a t /a a ) 2 (p + 2)(p + 4)(2+ph) 2 
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(p + 4)(ph + 2) R, 
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Here, h and p are the hollowing factor and index of the ini- 
tial density profile, ctj/cto is the initial space-charge inten- 
sity, and [£b / (2qX)](RiR' i )' is a parameter that measures 
the initial envelope mismatch of the beam. This nonlinear 
constraint equation can be solved numerically for fixed h, 
p, ai/<To and [£b/(2qX)](RiR' i )' to determine the ratio of 



final to initial rms radius of the beam (Rf/Ri). Employing 
the envelope equation (|j), the ratio of final to initial beam 
emittance is expressible as 

Rl 

Ri 



e xf 

C-xi 



I (Rf/Ri)* -[l-fa/ao)* 



(<7i/<j y 



R'l/{k%Ri 



"/30 J 



Eqs. (g|) and (^|) allow analysis of emittance growth from 
the thermalization of initial space-charge nonuniformities. 
We numerically solve Eqs. (|J) and (g) to plot (Fig. g) 
the growth in rms beam radius (Rf/Ri) and emittance 
(fixf/txi) due to the relaxation of an initial rms matched 
beam (R^ = = R") with nonuniform hollowed and 
peaked density profiles to a final uniform, matched pro- 
file. Final to initial beam ratios are shown for hollowing 
index of p = 2 and are plotted verses the "hollowing fac- 
tors" h (hollow initial density) and 1/h (peaked initial den- 
sity) for families of <Ji/(Jo ranging from di/ao — » to 
Ci/co —* 1- Growths are larger for the initially hollowed 
profile than the peaked profile and increase with stronger 
space-charge (smaller Oi/oo). However, the change in 
rms radius (Rf/Ri) is small in all cases, even for strong 
space-charge with strong hollowing (h — > 0) and peaking 
(1/h — > 0) parameters. Moreover, the increases in beam 
emittance (e x f/e X i) are surprisingly modest (factor of 2 and 
less) for intense beam parameters with Ci/ero ~ 0.1 and 
greater. At fixed (Ji/a^ and increasing steeping factor p, 
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Figure 2: Ratio of final to initial rms beam size (Rf/Ri) 
and emittance (e x f/e x i) verses h (a, hollowed beam) and 
1/h (b, peaked beam). 

similar modest growth factors are seen for hollowed beams 
for all but the most extreme hollowing factors (h ~ 0.1 and 
less), and as expected, much less growth is seen for peaked 
beams (closer to uniform). 

5 DISCUSSION 

The modest emittance growth at high beam intensities can 
be understood as general beyond the specific model em- 
ployed. Even for significant increases in emittance e x , the 
rms matched beam size is given to a good approximation 
by the envelope equation (Q) with the emittance term ne- 
glected. In this case, R ~ y/Q/k^o =const during the 
beam evolution and hence Rf ~ Ri. Employing the 
method of Lagrange multipliers, the free electrostatic en- 
ergy of the system at fixed rms radius (R) and line-charge 



(A) can be expressed as F = W — fd 2 x (/xir 2 + p-i)n 
with /ii.2 =const. Taking variations 5<f> subject to the Pois- 
son equation (0), one obtains to arbitrary order in 5(f), 



(9) SF = / d 2 x (q(/) - /iix 2 - fx 2 )Sn + 



9 |VM 2 



Thus, constrained extrema of F satisfy q<f> — pix 2 + /i2, 
corresponding to a uniform density beam centered on-axis. 
Variations about this extremum satisfy 5F > and are 
second order in 5<f>. Thus, the available electrostatic energy 
for thermalization induced emittance increase is modest for 
any smooth density profile. This can be demonstrated for 
our specific example using equation (Q) to plot AF — Wi — 
W f with Rt = R f verses h and 1/h for p = 2, 8 (Fig. g). 
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Figure 3: Free energy verses hollowing factors h and 1/h. 

It has been shown that the rms beam size and emittance 
undergo very small decreases on relaxation from a uniform 
density beam to thermal equilibrium over the full range of 
a./ov ( Min[e xf /e xl ] ~ 0.97 at a l /a ~ 0.45)0]. Thus if 
one views the relaxation as a multi-step procedure using the 
conservation constraints to connect the initial nonuniform 
profile to a uniform profile and then a thermal profile, any 
emittance growth will be experienced in the first step. This 
result together with the variational argument above show 
that the emittance growth results presented should be rela- 
tively insensitive to the form of the final distribution. 

Finally, caveats should be given for validity of the theory. 
First, the model assumes no generation of halo in the final 
state and that the initial nonuniform beam can be perfectly 
rms envelope matched. Initial mismatches can lead to halo 
production and provide a large source of free energy which, 
if thermalized, can lead to significant emittance growth[[l]]. 
Also, although the velocity space distribution is arbitrary in 
the present theory, choices that project onto broader spec- 
trums of modes will more rapidly phase mix and thermal- 
ize. Small applied nonlinear fields tend to enhance this re- 
laxation rate. Initial simulation results in a full AG lattice 
are consistent with model predictions presented here and 
will be presented in future work. 
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